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Oh" Karabulut and Sibert (J. Math. Phys. 38 (9), 4815 (1997)) have constructed 
<] 

\ an orthogonal set of functions from linear combinations of equally spaced Gaus- 

(N 

sians. In this paper we show that they are actually eigenfunctions of a q-oscillator 

O !■ in coordinate representation. We also reinterpret the coordinate representation ex- 



ample of q-oscillator given by Macfarlane as the functions orthogonal with respect 
to an unusual inner product definition. It is shown that the eigenfunctions in both 
q-oscillator examples are infinitely degenerate. 

I. INTRODUCTION 



> 

On 

cn 

f — ' Distributed Gaussians are a set of equally spaced Gaussians: g n (x) = e - c2 ( x - n ) 2 where 

O' 

(n = — oo, ...oo). A finite chain of them is often used in variational calculations as a basis 

•i— I , 

set. They are a very flexible and efficient basis set often yielding very accurate variational 
results. Calculating the potential matrix elements are often the most difficult part of a 
variational calculation. Because of their compactness it is very easy to calculate potential 
energy matrix elements with a few point Gauss-Hermite quadrature very accurately which 
is a major advantage of using this basis. 

We need orthogonal functions for variational calculations and the distributed Gaussians 
are not orthogonal. In a 1997 paper Karabulut and Sibert [ij] constructed a set of orthogonal 
functions from distributed Gaussians and they studied the underlying Gaussian quadrature. 
They called these functions distributed Gaussian polynomials (DG polynomials briefly). 
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Their motivation was entirely practical and they were not looking for an algebraic structure 
behind these function. Later Karabulut 2| showed how to normalize them and used these 
functions to construct a Wannier function set from distributed Gaussians. 

While searching for an operator that admits the DG polynomials as eigenfunctions the 
author came across a seminal paper by Macfarlane 3( in which he constructed a coordi- 
nate representation of his q-oscillator algebra, (the same q-oscillator was also studied by 
Biederharn J] but he did not have the explicit coordinate representation example that Mac- 
farlene gave). Eigenfunctions of Macfarlene were a linear combination of distributed Gaus- 
sians and it involved the q-binomial coefficients just like the DG polynomials. They looked 
similar to the DG polynomials but they were not the same. Following his example, the au- 
thor constructed another coordinate representation of the q-oscillator algebra that yields the 
DG polynomials as eigenfunctions. The q-oscillator turns out to be a coordinate represen- 
tation example of the Arik-Coon oscillator [5]. Macfarlane gave an inner product definition 
for his functions in term of Rogers-Szego polynomials. I also found a simpler inner product 
definition for his functions and reinterpreted his results. 

The outline of the paper is the following. In the second section I summarize basic results 
about the DG polynomials and discuss its links to Rogers-Szego polynomials. In the third 
section I derive the DG polynomials from the q-oscillator algebra. In the fourth section I 
discuss Macfarlane's example to his q-oscillator and reinterpret its eigenfunctions. Finally, 
in the last section I give a summary and discussion. 

II. DISTRIBUTED GAUSSIAN POLYNOMIALS AND THEIR PROPERTIES 

We will mostly denote the Gaussians e~° 2x2 as (f 1 where q = e~ c2 and c _1 is related to 
the width of the Gaussians. DG polynomials are defined as 

n 

$n(x) = Ec fc n (- 1 )V fc/ V x - fc)2 - (1) 

fc=0 

The CJ} are the well known q-binomial coefficients 

d = [ q ' q) \ , (2) 

{q,q)k{q,q)n-k 



where (q, q) n is defined as 

(q,q) n = (l-q)(l-q 2 )-(l-q n ), (3) 
and (q, q)o = 1. They satisfy the following orthogonality relation 2[]: 

§ n (x)$ m (x)dx = \\$n(x)\\ 2 5 nm , (4) 

where the norm ||$ n (x)|| is given as 



1/4 

2cV 

We will denote the normalized functions with lowercase 

3> n (x) a 



i^(*)ii = (^yv n/2 vW^ (5) 



£ C£(-l) Y^/Y*-*) , (6) 



where 



- 1/2 <2c^>* 



a 



We defined a this way for later convenience. 



Karabulut and Sibertfl] also found that the DG polynomials yield harmonic oscillator 
eigenfunctions in a particular limit as 

* n {s/y/2c) 2/2 

where h n (s) are the standard Hermite polynomials. We will refer to this limit later. 

DG polynomials are related to the Stieltjes-Wigert polynomials. Let us write the orthog- 
onality relation as follows 

oo 

/ $ n (x - s)$ m (x - s)dx = 0, (n^m). (9) 

— oo 

If we denote u = q~ 2x then <J> n (x — s) is written as 

$ n (x -s) = e~^ 2 ^- 4l ^u s P n (u; s), (10) 
where the polynomials P n {u; s) are 

p n ( U] s ) = jr c fc "(-i)y fe+s)2 - fc/ v. (ii) 

fc=0 
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Then the orthogonality relation becomes 

oo 

j e -(\nu?/(-2\n ri ) u 2s-lp n ^ U] s )P m ( u; s )& x = Q, ( n ^ m ). (12) 

— oo 

Evidently the polynomials -P n (w; s) are orthogonal with respect to the weight function 

W(U) = e -(ln«) 2 /(^ln,) u 2 S -l_ (13) 

For s = 1/2 the weight function is the lognormal distribution and the corresponding poly- 
nomials are known as the Stieltjes-Wigert polynomials. So the P n {u; 1/2) are proportional 
to Stieltjes-Wigert polynomials. 

The above connection to the Stieltjes Wigert polynomials were noted in Karabulut and 
Sibert 1]. Later Atakishiyev and Nagiyev 6j found a connection between the Stieltjes-Wigert 
polynomials and Rogers-Szego polynomials through the Fourier transform which implied 
that DG polynomials are also connected the Rogers-Szego polynomials. Here we note this 
connection. 

We take the convention for the Fourier transform as 



f(Q) = I ^ 6x f(x)dx. (14) 

— oo 

Using the Parseval relation of the Fourier transforms 

oo oo 

/ F*(x)G(x)dx= J F*(6)G(6)de j (15) 

— oo — oo 

we write the eq. @ as 

00 1/2 

J ^ n {e)<$> m {6)d6 = q- n (q,q) n 6 nm , (16) 



where $ n (#) is the Fourier transform of & n (x) given as 

$ n (0) = (1) V2 e -(^ 2 CJH-q-^e^y. (17) 



fc=0 



The polynomials 



H n {x) = Y,C n k x k (18) 

k=0 



are known as Rogers-Szego polynomials p, 13, |8|. Using them the orthogonality relation is 
written as 

00 / 2 \ V 2 

J H n (-q-V a e-™)H n {-q- 1 ^)e-W°r*M= f^J g)^- (19) 

—00 * ' 

A form of the Poisson summation formula reads 

J f(x)dx = J ( f; /(z + A;)W (20) 



v fe = — OO 

which is valid when J2 f(x+k) exists (in our case it does). Using this eq. ffT9l is expressed 
fc=— 00 

as 

.2 \ 1/2 



/ H n (- q -^e-^ e )H n (-q-^ e ) ( £ e" 

\k=—oo 



(21) 

The sum in parenthesis is periodic with period unity and it is a form of theta function. We 
can expand it in Fourier series. The Fourier coefficients can be calculated using the Poisson 
summation formula as 

00 / 2 00 

-2(7r/c) 2 (6>+fc) 2 _ / L n 2 /2 i27rn6» 



n 

Using the Jacobi $3 function |9(] defined as 



J2 e -2{n/cy{e+ky = r_ £ q n*/*cp*rd. (22) 

fc=— oo ' fc = — OO 



OO 

,3/ 



<W?) = E ? /2 e me , (23) 

n=— oo 

the orthogonality is written as follows: 
i 

J H^-q-^e-^H^-q- 1 ^ )^^ q)d6 = q- n (q, q)J nm . (24) 
o 

This relation is the well known orthogonality of the Rogers-Szego polynomials on the unit 
circle. Clearly it is the same thing as the orthogonality of the DG polynomials and one can 
be expressed in terms of the other. This is also of interest because Macfarlane 3j expressed 
orthogonality of his q-oscillator eigenf unctions in terms of orthogonality of Rogers Szego 
polynomials on the unit circle. 
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III. DG POLYNOMIALS AS A Q-OSCILLATOR EIGENFUNCTIONS 
A. Algebraic derivation of DG polynomials 



Lets define the translation operator T s as T s = e s sx . It has the effect of shifting a function 
to the left by s: T s f(x) = f(x + s). We define the creation and destruction operators a and 



at as 



at 



_L =T i/2 +i/4 _ T m (25) 
V 1 - Q 1 J 

' r g -+V4 _ T -l/2l T -l/2_ ( 2 g) 



Our inner product is the usual one 

oo 

U\9)= J r(p)g(x)dx, (27) 

— oo 

and the conjugate operator is defined as (f,Og) = (O^f, g). According to this (q x Y = q x 
and (^)t = — ^ and at given above is the right one. The a and at satisfy the commutation 
relation 

aat — gata = 1. (28) 

This is the commutation relation satisfied by the Arik-Coon oscillator Q]. 
We look for the eigenstates of the a? a 

a ] aA n {x) = X n A n (x). (29) 

We start from a 'ground state' A (x) that satisfies aA (x) = which yields the functional 
equation 

A (x + l/2) = q^ +1 ^A (x). (30) 

If we put Aq(x) = {cons.)w{x)q x2 in this equation we get w(x + 1/2) = w(x). Therefore the 
normalized ground state is 

2 

A (x) = a w w(x)q x , (31) 
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where w(x) is any (in general complex) function satisfying w(x + 1/2) = w(x) periodicity 
condition and CX •jjj IS the normalization coefficient 

-1/2 



ft,. 



|u>(x)| 2 q 2x2 dx 



(32) 



For w(x) = 1 we denote a w as just a whose value is given in eq. (J7]). We choose normalization 
of our eigenfunctions as (A n , A n ) = 1 and taking w(x) = 1 with this normalization will lead 
us to the normalized DG polynomials. 

Next we build the states (cft) n A (x) . Using the commutation relation one can easily show 
that if (a') n Ao(x) is an eigenfunction of a^a with the eigenvalue A n then the (a^) n+l Ao(x) 
is an eigenfunction of the a) a with eigenvalue A n+ i and one obtains a recursion relation for 
the eigenvalues 

A n+ i = q\ n + I. (33) 

Since aA (x) = 0, then A (x) is an eigenfunction of a) a with the eigenvalue Ao = and by 
induction it follows that all (a^) n A (x) are eigenfunctions. Using the recursion relation in 
eq. ( 1331) and Ao = we get the eigenvalues as 

1 -q n 



A, 



1-q 



(34) 



In exactly the same way that we do in solving the harmonic oscillator algebraically, we 
can easily obtain the following relations 

~ (35) 

(36) 



aA n (x) = y \ n A n _i{x), 
^A n (x) = ^X n+1 A n+1 (x). 



Then A„(x) can be written as 



A n {x) 



1 _ ri\ n 

s — ^(at)M (x). 



(37) 



Instead of applying n times, the following recursive relation is easier. Define A n (x) as 

a w w(x) 



A n (x) 



£D£(-l)ty n -*>/Y*- fc > a . 



(38) 



k=0 



Then if we apply a'/ y/X n +i to obtain A n+ i(x) and compare the coefficients we get the 
recursion relation for the D% coefficients. In this process w(x) completely commutes with 
because of periodicity: T^ 1 ^ 2 w{x) = w(x — 1/2) = w(x). The recursion relation we obtain 
is 

D n k +1 = q k D n k +D n k _ r (39) 

Together with the conditions D® = 1, D°_ 1 = this completely determines the D k . As 
can be shown easily, the q-binomial coefficients C k satisfy this recursion relation and the 
boundary conditions and therefore 

Dl = C n k = ^\ . (40) 

{q,q)k{q,q)n-k 

This shows that 

A^ x ) = a J W ^ jr C^-l) k q^' 2 q^\ (41) 

are the normalized eigenfunctions and they reduce to the normalized DG polynomials <p n (x) 
when w(x) = 1. 

Now let us see how the q-oscillator algebra reduces to harmonic oscillator algebra in the 
limit q — > (c — ■> 0). The limit in equation (jSJ shows us how to do it. First define the new 
variable z = ex. With this variable, the a will look like 



VI - e _u " 
and in the limit c — > this reduces to 



e (c/2)9/9z 



-c 2 /4 p -c Z _ p §(8/fte) 



(42) 



Similarly at reduces to 



1 9 

a (43) 



2*^-, (44) 



which, together with a, are the destruction and creation operators for the harmonic oscillator 
problem 

(-~ + Ai>(z) = Eil>{z). (45) 



Adz 



B. Further discussion on w(x) degeneracy 

Although the algebraic solution seems flawless, it is still very surprising that the A n (x) 
are orthogonal for all the functions w(x) satisfying w(x + 1/2) — w(x). Here we give a 
separate proof of it. 

When we multiply two of our parent Gaussians g n (x) = q( x ~ n ^ 2 we get daughter Gaussians 

G n (x) = g 2 ^-"/ 2 ) 2 as 

gn(x)g m (x) = q^ n - m)2/2 G n+m {x). (46) 

The parent Gaussians are centered at integers whereas the Daughter Gaussians are centered 
at both integers and half integers. Therefore the (f) n (x)(f) m (x) product of the normalized DG 
polynomials 

n 

4> n (x) = -j^= £ CU-l) k q in - k)/2 q {x ~ k) \ (47) 
^{q,q)n k=o 

can be written as a linear combination of Daughter Gaussians as 

n+m 

<t>* n {x)ct> m {x) =a 2 J2 drq 2{x ~ k/2)2 ■ (48) 

fc=0 

If we integrate this we should get 5 nm due to the orthogonality of normalized DG polynomials. 
The integrals Jf^, q 2 ( x ~ k ^ 2 dx are independent of k (just shift the integral by k/2) and have 
the value 1/a 2 . Therefore we get the relation 

n+m 

£ dr = Sum- (49) 

k=0 

Now, the A* n (x)A m (x) product of the A n (x) functions 

Aji ( x ) = a J Wi - X) CZ(-l) k q {n - k)/2 q (x - k) \ (50) 
^J(q,q)n k=o 

can be expressed as 

n+m 

A* n (x)A m (x) = \a w \ 2 J2 d n k m (\w(x)\ 2 q^- k ^ 2 ) . (51) 
If we integrate this we get 

J A* n (x)A m (x)dx = £ d n k m [ \a w \ 2 J \w{x)\ 2 q 2{x - k ^ 2 dx \ . (52) 



k=0 



k=0 



-oc 
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The integral in parenthesis can be shifted by k/2 as 

oo oo 

J \w{x)\ 2 q 2{x - k/2)2 dx= J \w{x + k/2)\ 2 q 2x2 dx. (53) 

— oo — oo 

Because of the periodicity of w(x), we have w(x + k/2) = w(x) and all the integrals are 
independent of k and they have the value 1/ \a w \ 2 . This yields 

J A* n (x)A m (x)dx = ]T dr = 5 nm} (54) 

— oo k=Q 

which follows from eq. (l4"9]) . Therefore the orthogonality holds for any w(x) satisfying the 
periodicity requirement. 

We actually found more than the DG polynomials from the algebraic treatment. We 
found and infinite set of orthogonal functions. 



IV. A DIFFERENT INTERPRETATION OF Q-OSCILLATOR EXAMPLE OF 

MACFARLANE 

In a seminal and widely cited paper Macfarlane constructed a different coordinate repre- 
sentation of the q-oscillator. His definition of the creation and destruction operators are 

b = e 2x - e x e s{d/dx \ (55) 

Obviously, according to the usual definition of the inner product in eq. (1271) the b^ is not 
right. Macfarlane himself goes on to construct eigenstates of tfb without discussing the inner 
product and orthogonality relation first. After finding the eigenf unctions he states inner 
product and orthogonality relations in terms of orthogonality of Rogers-Szego polynomials 
on the unit circle. The functions we find do not agree with the functions he found and 
apparently his formula for the eigenf unctions contains some error. Therefore we redo the 
problem with a different inner product definition here. 
We define the inner product as 
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(f,g) = J (Pf*(x))g(x)dx. (57) 

— oo 
oo 

= / f*(-x)g(x)dx (58) 

— oo 

where P is the usual parity operator. For an operator O, its conjugate is defined by the 
relation 

(f,Og) = (&f,g). (59) 

The &t expression in eq. (1561) is right if one defines the inner product this way because 
(e x y = e~ x and (4;)^ = 4:- Notice that this definition of the inner product does not 
satisfy (/, /) > condition. But we will not use this property of the inner product in our 
development. We will use the conjugacy definition in eq. (!59|) in a few places. 

In order to obtain functions as a linear combination of the Gaussians centered at non- 
negative integers we change the variable x = —c 2 (y + 1/4) and s = — c 2 /2 and we take 
q = e~° 2 as before. We will also divide b and by \Jq{l — q) which does not change the 
eigenfunctions (only eigenvalues), but in the limit of the harmonic oscillator (q — > 1) it helps 
to get the right results. The new b and W are 

b = (q 2y+1/2 - g^/4 T i/2) /^(TT^ (60) 

P = (V _ T l/ V , + l/4) /JfrZjj, (61) 

where T 1//2 = q^^v)/ 2 . They satisfy the commutation relation 

b%-qW = l. (62) 

Note that this is somewhat different than the commutation relation in eq. (1281) . 
Again we want to find the eigenfunctions of the operator tfb 

VbB n {y) = \ n B n {y). (63) 

We again start from a ground state that satisfies bB (y) = (hence A = 0). bB (y) = 
yields the same functional equation that Aq(x) satisfy 

B (y+ 1/2) =qy +1 / 4 B (y), (64) 
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which we already know has the solution B {y) = a w w{y)q y where w(y) satisfies the w(y + 
1/2) = w(y) periodicity condition. The infinite degeneracy of the states appear here too. 

Next we build the unnormalized states (6^) n i?o(l/)- Using the commutation relation we 
can easily show that if the (6^) n i?o is an eigenfunction of tfb with an eigenvalue A n then 
(6^)™ +1 i?o is an eigenfunction with the eigenvalue A n+ i where A n+ i is related to the A n as 

q\ n+1 = A n - 1. (65) 

Since B (y) is an eigenfunction, by induction all (tf) n B are eigenfunctions too. Starting 
from Ao = 0, this recursion relation yields 

A„ = -q~ n {j^j • (66) 

Notice that the eigenvalues are negative. 

To obtain relations similar to eqs. (l35|l36l) We form the inner product 

\ n (B n , B n ) = (B n ,PbB n ) = (bB n ,bB n ). (67) 

Taking bB n = \i n B n _\ we get 

■v I i2 (B n -i, B n _i) 

Xn = 1^1 (n n \ ■ ( 68 ) 

Since |/i n | 2 is positive and A n is negative, the (B n _i, B n _i) and (B n , B n ) must have opposite 
signs. This looks surprising but we should remember that in our definition of the inner 
product the (/, /) can be negative. Therefore we will take our normalization as 

(B n ,B n ) = (-l) n . (69) 

With this normalization we obtain /i n = y/—X n : 



bB n = sj-KB^. (70) 

Starting from 



An+i-Bn+i = WbB n+ i = J—X n+ i¥B n , (71) 
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we also obtain 



b^B n = —J—\ n+ iB n+ i. (72) 



This result is also a little unusual because of the sign in front. It is a consequence of negative 
A n eigenvalues. 

Since we started from an unusual inner product definition there might be doubts on 
orthogonality of the B n . Consider the inner product 

(B m ,PbB n ) = (PbB m ,B n ). (73) 

Using b^bB k = X k B k we get 

(X n -X m )(B m ,B n )=0, (74) 

which yields (B m , B n ) = when m ^ n. We added this common proof to emphasize that 
the orthogonality does not depend on the (/, /) > property of the usual inner products. 
We just used the conjugacy relation eq. (l59l) in this proof of orthogonality. But unlike the 
case of usual definition inner product given in eq. (l27l) . (/, /) = does not imply / = when 
we don't have (/, /) > property. 

To complete the discussion we obtain the eigenfunctions B n (y). It makes things easier to 
figure out the coefficient of Gaussian centered at zero(q y2 ) first. The 



V = (g- 2 ^ 2 - rVYiH-V*) / v^O^), (75) 

has two parts that do different things. When g~ 2 f +1 / 2 acts on g(^~ fc ) 2 it produces the next 
Gaussian g(f _fe_1 ) 2 and multiplies it with a constant. When T x l 2 q~ yJrX l A act on q^ y ~ k ^ 2 it 
gives q( y ~ k } 2 back and multiplies with a constant. Here are the precise relations 

q -2y+l/2 q (y-k) 2 = q -2k-l/2 q (y-k-lf ^ ^ 
T l/2 q -y+l/A q (y-k) 2 = q -k q ( y -k) 2 _ ^ 

From the second relation we have T 

l/2 q -y+l/4 q y 2 = g y 2 

which means this operator leaves 
q y2 as it is. The other operator g~ 2 ^ +1 / 2 creates k = 1 Gaussian (g^ -1 ^ 2 ) from it. Since we 
produce higher eigenfunctions by applying —tf / y/—\ n successively, after each application 
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the coefficient of k = Gaussian changes by a factor l/yg(l — g)(— A n ). We denote the 
coefficient of k = Gaussian in B n (y) by and it should be 

Cn C^UI / ; = &W ; = • ( 78 ) 

vV(i-<?r ^/(-Ai)(-A 2 )...(-A n ) yW)^ 

Now let us take the B n (y) of the form 

n 

B n (y) = w(x)( n Y / Ett y - k) \ (79) 

k=0 

where Eq = 1 by construction. We can use Eq. ( J72|) to generate a recursion relation for E% 
as before. To show a different and easier way of doing things we will use eq. (1701) relation 
this time. The b has the effect of shifting each Gaussian to the left by one unit 

l q (y-k? = _ g (1-9 J gfr-fc+i)^ (g0) 

v^l - ?) 

and it also destroys the leftmost Gaussian (the = Gaussian centered at zero). Therefore 
the equality 



b \{CnE n k ) q^ 2 ] = ^K. [(Cn-lK-1 1 ) g fe " fc+1)2 ] (81) 

should hold for each Gaussian. This yields the recursion relation for the E% as 

El = -Et_l (j^j q~ n ~ k+3/2 . (82) 

Together with the condition Eq = 1 this is enough information to solve the E%. The 
quotient in parenthesis tells us that the q binomial coefficients C£ are involved. If we set 
E% = (— l) k C%q u ( n ' k ^ we get the recursion relation for u(n, k) as 

u{n, k) - u(n - 1, k - 1) = -n - k + 3/2. (83) 

Together with the condition u(n, 0) = (which follows from Eq = 1) this is uniquely solved 
as u(n, k) = —nk + k/2. Therefore the B n (y) should be 

„n(n— 1)/4 n 

B n (y) = a w w(x)^= £ C£(-l) V^N^- (84) 
\J{q,q)n k=o 
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This formula does not agree with the result of Macfarlane even after setting c 2 = —2s 
and y = —l/4 + x/2s back in B n (y) above. To be sure that we have the right formula we 
have checked numerically if the B n (y) satisfy the orthogonality relation 

J B n (-y)B m (y)dy = (-l) n 6 nm . (85) 

— oo 

They satisfy it perfectly and we are sure that we have the right formula. Apparently Mac- 
farlane's paper contains an error. 

The harmonic oscillator limit is as straightforward as it is in DG polynomials case. We 
change variable y = z/c and take the limit c — > which yields 

S ^*-53S- (86) 
p - ->- \l < 87 > 

Then WbB n = X n B n reduces to (A n — > —n in this limit) 

(-\^ + z2 )*n = {n + l/2)^ n . (88) 

The fact that the inner product is defined differently makes no difference in this limit because 
the harmonic oscillator eigenfunctions are either even or odd. The harmonic oscillator 
eigenfunctions satisfy the normalization condition 

J ^ n (-y)^ m (y)dy = (-l) n 5 nm (89) 

as can be verified easily using parity of the wave functions. 

Now, just as the DG polynomials, the orthogonality of the B n (y) can be expressed as 
an orthogonality relation of the Rogers-Szego polynomials on the unit circle. We will take 
w(x) = 1 and a w = a for this. By Fourier transforming the orthogonality relation eq. (|85p 
we get the relation 

oo 

J B* n (-9)B m (6)d6 = (-l) n 5 nm , (90) 

— oo 

where B n (6) is the Fourier transform of the B n (y): 

/ -7T \ 1/2 n 

B n (9) = ( - 2 ) ( n e-^ 2e2 £ CU-q- (n - 1/2) e^y. (91) 



fc=0 
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Here 6 is the Fourier transform variable just as before. Using the Poisson summation formula 
in eq. (120]) the orthogonality relation can be expressed as 

/ H n (-q^-V^ e )H m (-q-^-V%™) ( £ e -W*r<fi+>A d6 = - { -^5 nm . (92) 

Setting the sum in parenthesis from the eq. ff22|) and fl23|) we obtain 
i 

J ^ n (-g-("- 1 / 2 )e i2 ^)iJ m (-g-( m - 1 / 2 )e i2 ^)^(2^;g)d0= ^^^(g, q) n {-l) n 8 nm . (93) 
o 

This is again some form of the orthogonality relation of the Rogers-Szego polynomials on 
the circle. This is a new set of orthogonality relations and we are not aware of its existence 
in mathematical literature. 



V. SUMMARY 

In this study we found that the DG polynomials that Karabulut and Sibert discovered 
before are actually eigenf unctions of coordinate representation of the Arik-Coon q-oscillator. 
We derived the DG polynomials from q-oscillator algebra. We also indicated that orthogo- 
nality of the DG polynomials can be cast into orthogonality of the Rogers-Szego polynomials 
on the unit circle and the two are equivalent. 

We showed that the example given by Macfarlane can be interpreted with an unusual 
inner product definition and we constructed the corresponding orthogonal functions. Their 
orthogonality can be recast into a form of orthogonality relation for the Rogers Szego poly- 
nomials on the unit circle. We were not able to find this result in mathematical literature 
and it is probably a new result. 

A very interesting result of this work is that the eigenstates of the q-oscillators we solved 
turned out to be infinitely degenerate. We know that algebraic solution of the one dimen- 
sional harmonic oscillator is nondegenerate. In the usual algebraic solution of the harmonic 
oscillator we base our arguments on the commutation relations and the algebra has nothing 
in it that implies nondegenerate states. But the ground state 20 o = yields a unique solu- 
tion for the harmonic oscillator because it is a differential equation and nondegeneracy of the 
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excited states follows from this. For our q-oscillator examples we have a first order difference 
equation for the a<f)Q = and such equations together with the boundary conditions do not 
uniquely define a function. It defines a function on all real axis if its values in a 1/2 wide 
interval are known. 

Finally, the freedom to choose w(x) arbitrarily gives us possibility of constructing or- 
thogonal function sets more general than the DG polynomials. Consider a set of func- 
tions Wo(x), Wi(x), w 2 (x).... satisfying the periodicity condition w n (x + 1/2) = w n (x) 
(n = 0, 1, 2, ...) and the orthogonality relation 

CO 

J w* n (x)w m (x)q 2x2 dx = 5 nm . (94) 

— oo 

Then the set of functions 

r„ m (x) = a ( Wn ) Wn (x)(f) m (x) (95) 
a 

will satisfy an orthogonality relation of the form 

CO 

J T* nm (x)T ij (x)dx = 5 ni 5 mj . (96) 

— oo 

Here a{w n ) is the a w for w = w n (x) given in eq. (|32l and a is given in eq.([7j). This gives 
us much freedom to construct orthogonal function sets useful as basis sets in variational 
calculations. This possibility should be investigated in future research. 
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